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SUMMARY 

Meteorological and dynamical requirements pertaining to the specification 
of middle atmosphere tides by the MST radar technique are outlined. Major is- 
sues addressed include; ( 1 ) the extraction of tidal information from measure- 
ments covering a fraction of a day; ( 2 ) the ramifications of transient effects 
("tidal variability") on the determination and interpretation of tides \ (3) re- 
quired temporal and spatial resolutions and; (4) desired global distributions of 
of MST radars, so as to complement existing MST, meteor wind, and partial re- 
flection drift radar locations. 

Specific conclusions are as follows: 

(1) A height resolution of 2 km is considered necessary to resolve all tidal 
structures that might reasonably be expected throughout the middle atmosphere 
( 10-100 km). 

(2) Time resolutions of about 30 minutes in the stratosphere and 5 minutes in 
the meso- thermosphere are required to determine with acceptable accuracy tidal 
Fourier components from a single day's measurements. (It is assumed that data 
above 50 km are available only during daytime for 12 hours.) Alternatively, one 
to three weeks of data at 30 minute intervals would enable determination of 
semidiurnal and diurnal components above 50 km to roughly the same accuracy. 
(These estimates assume nominal values of diurnal and semidiurnal amplitudes and 
standard deviations relative to the standard deviation of uncorrelated noise in 
the data, and may require adjustment to reflect actual conditions.) 

(3) Nominal accuracies necessary for the delineation of middle atmosphere tidal 
structures are about 2 m sec~l (.5 m sec - ^) for altitudes above (below) 60 km. 

(4) Since the set-up time for tides is typically several days or more, it is 
not meaningful to speak of "day-to-day variability" of a tide (which implies a 
steady-state response at a period exactly a subharmonic of a solar or lunar 
day). Four-day vector averages are recommended for comparison with models, 
other data, and delineating the variability of tidal components. 

(5) The geographical distribution of radars is an extremely important factor in 
determining the scientific return of global tidal measurement. 

( 6 ) Data base management and standardization of reduction and analysis 
techniques are extremely important issues with regard to tidal studies, which 
require coordination and interpretation on a global scale. 

INTRODUCTION 

Atmospheric tides, oscillations in meteorological fields at subharmonics of 
a solar and lunar day, comprise an important component of middle atmosphere 
dynamics. Solar tidal components observable by MST radars are primarily excited 
by insolation absorption by H 2 0 and O 3 in the tropospheric and strato- 
me 8 ospheric regions, respectively. Tidal oscillations in the middle atmosphere 
can to first order be viewed as the superposition of several "quasi-modes" each 
with somewhat distinguishable and identifiable horizontal structures and verti- 
cal wavelengths. The modifier "quasi" is used here to suggest possible devia- 
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tions from "true” model structures (that is, corresponding to the eigenfunctons 
of Laplace's tidal equation) due to transient effects or the influence of lati- 
tudinal ly varying background propagation conditions. The specific mixture of 
these structures, which might be comprised of migrating and non-migrating, ver- 
tically propagating and evanescent, and symmetric and asymmetric components, de- 
termines the global spatial variability of the tidal oscillations. 

The purpose of this report is to outline meteorological and dynamical re- 
quirements pertaining to the specification of upper atmosphere tides by the MST 
radar technique. In a report (LINDZEN , 1969) prepared for the Panel on Dynamics 
and Structure of the Neutral Atmosphere of COSPAR's Working Group II, a group of 
scientists addressed the problems of measuring tides and gravity waves in the 
upper atmosphere. HAURWITZ (1964) is referenced extensively in their report. 
Some of their main points are as follows: 


(1) A minimum of five daily rocket ascents, preferably equally spaced in time, 
are needed to determine diurnal and semidiurnal tidal components at a given 
station (HAURWITZ, 1964). 

(2) Statistical measures of the reliability of tidal measurements, for instance 
as given by the radius (R) of the probable error circle, should always be pro- 
vided (HAURWITZ, 1964). 

(3) Quoting HAURWITZ (1964), "According to error theory R decreases as 1/ /N 
where N is the number of determinations. In order to reduce R to half its value 
one needs four times as many determinations, that is a four times longer series. 
In this manner it will be possible to estimate even from a fairly short series, 
having a large R, how much longer the observation program should be continued". 

(4) The question of how many rockets must be deployed in order to measure 
tides is difficult to answer. It involves a tradeoff between cost and the 
degree of reliability one is aiming for. 

(5) Tides can be associated with rather short vertical wavelengths, as low as 
12 km for important diurnal components, A vertical resolution on the order of 
at least 3 km is needed so as not to smooth out such components. 

(6) Meteor radars appear capable of routinely providing economical measurements 
of tides with good spatial and temporal resolution. 

(7) Observations may never be sufficiently extensive to properly define upper 
atmosphere tidal oscillations. A hybrid approach involving both theory and 
experiment may be necessary to arrive at a comprehensive delineation of the 
phenomenon. 

(8) Valuable consistency checks on theory can be performed with less data than 
would be necessary for a complete tidal determination. 

The above report emphasizes the role of rocket soundings in detecting upper 
atmosphere tides, but their points remain valid in the context of radar measure- 
ments, Note that they did recognize the potential role that meteor radars could 
play in the observational study of tides. A shortcoming of the meteor radar 
technique is that it is restricted to the 80-100 km height range. This defi- 
ciency is offset by the global array of meteor radars and partial reflection 
drift radars that now exist for dynamical studies. The MST radar can in princi- 
ple provide wind measurements over the 10-100 km height range, and thus add an- 
other dimension to the experimental determination of tides. 


There have also been significant advances in state of the art modelling as 
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well as upper atmosphere measurements techniques since the report by Lindzen and 
his colleagues. In this light the present report defines operational require- 
ments of MST radars necessary to provide useful tidal information, and assesses 
the overall role of the MST radar technique in the future of tidal studies. 

MEASUREMENTS COVERING A FRACTION OF A DAY 

In the mesosphere and lower thermosphere nighttime ionization levels are 
often so low as to preclude wind measurements by the MST or Thomson scatter ra- 
dar techniques. The question then arises as to the exact limitations imposed by 
deriving tidal components from data covering only a fraction of a day. Follow- 
ing a recent treatment by CRARY and FORBES (1983), assume that the data consists 
of mean, diurnal, and semidiurnal components with a random noise component of 
standard deviation a. Consider a fit function of the form: 

y( t ) = b Q + b ± sin t + + ) (1) 


and define to be the standard deviations of the fit coefficients. In the 
case of a full-period (24-hour) fit of the above function it is well known (see, 

2 

for example, SCHUSTER, 1898) that is given by « a b2 ” Tjjp, where N is 

the number of points fit. Thus, for this case, a fit to even a small number of 
data points (i>;5) reduces the error in the fit coefficients to below that of the 
data points. In the case of a partial-period (< 24-hour) fit, however, the 
errors in the fit coefficients can be large enough, even for relatively small 
fluctuations, to make the extraction of meaningful results from the data 
unfeasible. This problem is briefly discussed by CHAPMAN and LINDZEN (1970). 

In order to illustrate the preceding statements and give a rough estimate of the 
errors involved in an actual fit, the coefficients a bi were obtained numerically 
from random data with o » 1. It is simple to show that the dependence on a is 
linear: 

a M = f. (T,N)o i * 0,1,2 (2) 

where T * fit span and N - number of points fit. Figure I shows f^T/N) for 
N = T + 1 (hourly data) for 1*7 through T ® 20. For large values of T the 
fj(T,;N) show the expected convergence fi(24, 25) « .4 for i = 1,2. For smaller 

values of T, however, the various f,(T,N) show a rapid increase with decreasing 

T; fjXT/N) * 1 (a bi * a) for T ^ 15 hours. For T - 12 hours a^ * 3.1a and 
a b2 * Figure 1 also shows f^(T,N) for N = 4T + 1 (data quarter-hourly). 

In this case ^(24, 97) = .2 and f^TjN) ~ 1 for T * 12-14 hours. 

Table 1 contains data for T « 12 hours (a half-day’s measurements) for 
various equally spaced time intervals t between measurements (t (minutes) = 

6QT/ (N-l) ). Be aware that a^ or a t 2 * o is not necessarily a good absolute 
criterion for acceptability, but is only a relative measure. For a given a 
reflected in the data, one must choose N and/or T such that the a^'s are 
acceptable. Let us take a typical example. Suppose we have 12 hours of daytime 
wind measurements at a given altitude. We bin these data in local time with a 
time interval (t) that is as yet unspecified. A Fourier fit to the data is 
performed and there exists a noise component in the data characterized by some 
standard deviation, a. For the sake of argument, assume b^ * 10a and b 2 85 5a 
and that we require ^bl/bl and a b2^b2 to be £* 1# 90 ^ aa altitude a * 5 m 

sec” 1 might be a typical value.) It is evident from Table 1 that these 
constraints require t $ 5 minutes. Similarly, for b^ * 5a and b 2 = 10a we 
obtain t < 2 minutes to measure both the diurnal and semidiurnal components to 
the desired accuracy. To satisfy only the semidiurnal constraint in this latter 
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Figure 1* vs T for hourly data (N*T+1) 

and quarter-hourly data (N=4T+1) where N = 
number of data points and T = fit span 
(hours) . 


example would require t < 20 minutes. 

Alternatively, we can examine how many days (N) of data for a larger x are 
required to obtain the same accuracies. For x =* 30 minutes, and the case where 
b^ * 10a and b 2 * 5 a, a * 2.19a and a^ * 1.27a for one day's measurement. 

To obtain a^, 1 a we require 2.19/vfti “ 1> and for 0 b2 * *5a'we require 
1.27//N2 ** *5; therefore Hi « 4,8 and N2 88 6.5. To satisfy both constraints 

we must take the larger of and N 2 , or N * 7. Similarly, when b^ = 5a and 
and b 2 “ 10 a , 2.19/ *4?-. * 15 and 1.27/ No * 1; hence Nn = 19.2 and N 2 * 1*6, 
so that N * 20. 


Table 1. Standard deviations of Fourier components 
for 12 -hour fit span 



Data 

Interval (x 

, minutes 

> 



(60) 

(30) 

(15) 

(5) 

(2) 

c bl /0 

3.10 

2.19 

1.55 

.89 

.57 

ro 

Q 

1.80 

1.27 

.90 

.52 

,33 
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It must be emphasized that these results are derived for noise which is 
uncorrelated over the time interval between data points. Actual data, of 
course, may or may not reflect this; the above is intended only as a guideline 
for minimum criterion for reliability. 

VARIABILITY 

Observations of atmospheric tides are often characterized as being variable 
on a day-to-day basis. However, BERNARD (1981) points out that if the intrinsic 
time (t s ) to set up stationarity around the earth for a particular mode is 
greater than a characteristic time scale associated with the variability, then 
transient effects will exert some influence over the observed 1 tidal* behavior. 
Therefore, it may only be meaningful to refer to the variability of a tide over 
time scales greater than x s . Furthermore, theoretical models should be compared 
with a steady state or "average” tidal component uncontaminated by transient ef- 
fects. 


As a guide to interpreting observations, BERNARD (1981) has calculated 
characteristic propagation times (x p ) around the earth of various pertinent 
tidal modes (see Table 2). These are obtained from horizontal group velocities 
(VgH, a function of latitude) determined by an "equivalent gravity wave (EGW)" 
approximation (SPIZZICHINO, 1969; RICHMOND, 1975) wherein the value of Vg H is 
evaluated the latitude (0 O ) where the vertical group velocities of the tidal 
mode and EGW are equal. These propagation times provide some measure, at least 
a relative one, of the set-up times (t s ) for various modes. Note that x p is 
on the order of ten days for the (2,6) and (1,1) modes, and hence the condition 
of stationarity may not be frequently realized in practice. 

It is difficult to draw specific quantitative information from these 
estimates, A recent unpublished calculation by Forbes indicates that for 
x y - 3-5 days, distortions of the (1,1) modal shape are typically 15% in 
amplitude and 1 hour in phase, suggesting that the tabulated values of x p may 
be too high by more than a factor of two with regard to estimating the onset of 
"measurable" effects. Given the uncertainties inherent in these estimates, it 
appears reasonable to say it is meaningful to refer to variations in most tides 
over periods of over 4 days, whereas variations over shorter time scales may 
reflect transient effects. To achieve the highest degree of consistency between 
data sets and provide a meaningful comparison with theoretical models, it is 
therefore recommended that tidal measurements from MST radars be represented as 
vector averages of no less than 4 day* s observations. 


Table 2 .* 


Mode 

0 Q (deg) 

V (km hr ^) at 0 
gh o 

T p (days) 

(2,2) 

35.4 

662.4 

2.50 

(2,4) 

43.3 

273 . 6 

6.08 

(2,6) 

46.4 

165.6 

10.1 

(1,1) 

21.2 

140.4 

11.9 

(1,3) 

21.1 

61.2 

27.3 


*After BERNARD (1981), See text for definitions of 9 , V t . 

o gh p 
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VERTICAL RESOLUTION AND ACCURACY 

Vertical wavelengths associated with tidal modes expected to be observable 
in the mesosphere and lower thermosphere (50-100 km) generally fall in the range 
30-70 km. (See FORBES, 1983, for a list of references.) Values as low as 20 km 
associated with diurnal propagating tides might be found at tropical latitudes 
( <30°) . However, there exists evidence from observations and theoretical models 
that the superposition of modes, particularly in the 80-100 km region, can 
cancel and reinforce so as to produce structural features with scales of the 
order of 10 km. 

In the stratosphere, observations (GROVES 1980; FUKAO et al., 1978; and 
WALLACE and TADD, 1974) are frequently characterized by short vertical wave- 
lengths (~ 10-20 km), particularly in connection with the diurnal component. 

KATO et al. (1982) investigated the generation and propagation of tides due to 
geographically localized sources of excitation, and are able to explain many 
characteristics of the observed structures, which are apparently non-migrating 
tidal components. Stratospheric tidal amplitudes are typically much smaller 
than in the mesosphere, and may reflect transient effects due to variations in 
source characteristics. Both of these factors suggest that data series of a 
week or more may be necessary to properly define the stratospheric tidal compon- 
ent, assuming time resolutions of 15-30 minutes and 24 hours of measurements per 
day. 


The above information indicates that a vertical resolution of 2 km through- 
out the 0-100 km height region would be necessary to resolve all tidal struc- 
tures that might reasonably be expected. In addition, observations suggest that 
nominal accuracies necessary for the delineation of middle atmosphere tidal 
structures are about 2 m sec"^ (.5 m sec - -^) for altitudes above (below) 60 km. 


RADAR NETWORKS 


Our understanding of middle atmosphere tides can best be improved and 
tested with simultaneous global measurements, preferably covering a week or 
more at a time. It is important to span a range of latitudes within a longitude 
sector, or vice-versa, in order that modes (migrating and non-migrating) be sep- 
arated. Insofar as it is feasible, the placement of newly proposed MST radars 
should complement existing MST, meteor wind, and partial reflection drift radar 
locations. Tidal theory and theoretical simulation models of tides can play an 
important role in defining the locations and operating requirements of radars so 
as to derive the optimum scientific return. 

Ideal tidal campaigns would include observatories geographically conjugate 
about the equator and at least one equatorial radar. A MAP Study Group under 
the direction of Prof. S. Kato is currently examining the scientific aspects of 
an equatorial radar. Recently, geographically symmetric radars at Kyoto (35 °N) 
and Adelaide (35 °S) have delineated strongly asymmetric tidal behavior about the 
equator in both diurnal and semidiurnal components, between 80 and 100 km 
altitude (ASO and VINCENT, 1982). For "true" modes (that is, eigenfunctions of 
Laplace's tidal equation), asymmetric tides are characterized by nonzero N-S 
velocities and zero E-W velocities at the equator, whereas the reverse is true 
for symmetric components. A properly placed equatorial radar could, therefore, 
contribute significantly to the deconvolution of tidal structures in the middle 
atmosphere. 

Tidal observations at the equator could also contribute significantly to 
other aspects of middle atmosphere dynamics. For instance the gravest symmetric 
diurnal propagating (1,1) mode is thought to be capable of producing significant 
turbulence via nonlinear cascade or convective instability mechanisms, and a 
mean easterly jet (~50 msec"^) peaking near 105 km due to momentum deposition 
in the region where the tide is dissipated (MIYAHARA, 1981). An equatorial 
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radar would, therefore, also provide the opportunity to observe the important 
middle atmosphere processes of wave-mean flow interaction and wave/turbulence 
coupling for a known wave forcing . However, these phenomena primarily occur 
above 75 km, so that it is imperative that measurements extend as far into the 
lower thermosphere as possible, preferably above 105 km. Therefore, placement 
of the radar must avoid possible interference by the equatorial electro jet. 

Given the necessity of global coordination, it is important that some 
consensus be reached concerning the appropriate method of extracting tidal 
components from data time series, expressing an appropriate Average* tidal 
structure for comparison with models, and describing variability about mean 
monthly or mean seasonal behaviors. In addition, the importance of data ex- 
change and data base management cannot be overemphasized for the global study of 
tidal phenomena. 
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